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The envelope-function method with generalized boundary conditions is applied to the description 
of localized and resonant interface states. A complete set of phenomenological conditions which 
restrict the form of connection rules for envelope functions is derived using the Hermiticity and 
symmetry requirements. Empirical coefficients in the connection rules play role of material param- 
eters which characterize an internal structure of every particular heterointerface. As an illustration 
we present the derivation of the most general connection rules for the one-band effective mass and 
4-band Kane models. The conditions for the existence of Tamm-like localized interface states are 
established. It is shown that a nontrivial form of the connection rules can also result in the forma- 
tion of resonant states. The most transparent manifestation of such states is the resonant tunneling 
through a single-barrier heterostructure. 



I. INTRODUCTION 

Over 50 years the effective mass or envelope function 
method is widely used to describe physical properties 
of various spatially inhomogcncous semiconductor sys- 
tems. Originally the envelope function approach was 
developed for external potentials which vary slowly on 
the atomic scale 1 . Nonetheless the application of this 
method to semiconductor nanostructures with micro- 
scopically abrupt heterointerfaces is commonly accepted 
and frequently gives unexpectedly good results. 2 

The central question to the effective mass method in 
the context of heterostructure applications is how to con- 
nect envelopes at a heterointerface. The simplest and 
historically the first connection rules follow from the as- 
sumption that the effective Schrodinger equation for the 
envelope function ^(r) 

Hip(r) = £ty(r) (1) 

is valid everywhere in space, and any abrupt variation of 
material parameters can be viewed as a limit of smooth 
function. This assumption allows to obtain boundary 
conditions by the integration Eq. (1) over an infinitesi- 
mally small distance across the interface 3 . As a result we 
get the connection rules 2,4 

V(+o) = v(-o), 

OH 

vM+o) = M»(-o), fi« = ^— , (2) 

(p z is the component of momentum perpendicular to the 
interface) which are commonly called the standard or 
BenDaniel-Duke boundary conditions. For the one-band 
effective mass model the boundary conditions Eq. (2) are 
reduced to the continuity of an d —d z ip, where m is the 
effective mass. Obviously, the standard connection rules 
cannot be universal since they contain only bulk parame- 
ters of materials which constitute the heterojunction and 
thus completely neglect internal properties of the inter- 
face. In fact, they work quite well for GaAs/AlGaAs het- 
eropair, but fail to describe properly more complicated 



situations. 5 ' 6 The failure of the standard boundary con- 
ditions was also explicitly demonstrated for a number of 
particular microscopic models of heterostructures. 7 ' 8 ' 9 ' 10 
A natural phenomenological way to take into account 
the above mentioned results is to relax the assumption of 
applicability of Eq. (1) near an interface and to allow for 
the discontinuity of both envelops and their first deriva- 
tives. Due to the superposition principle, wave functions 
at opposite sides of the interface must be connected by 
a linear relation. Assuming locality of this relation we 
arrive at the following connection rules 

( ^(+0) \_(Tu T 12 \ ( V(-0) \ 

^ dM+o) )~ V T 2i T 22 J \ 0»V(-o) J ( ' 

where = T is the transfer matrix. Boundary con- 
ditions of the form Eq. (3) were introduced by Ando 
and Mori for the one-band effective mass model. 5 Later 
they were also adopted to different milti-band models 
which describe V and X conduction bands, 11 ' 12 degen- 
erate heavy hole and light hole valence bands, 13 and 
conduction and valence bands within the spherical Kane 
approximation. 1 

The transfer matrix approach unifies all possible 
boundary conditions which have been suggested in the 
literature (for particular examples see Refs. 7,15,16,17, 
18,19). In fact, the standard boundary conditions Eq. (2) 
are described by the diagonal transfer matrix: T\\ = 1, 
Tii = T 2 i = and T 22 = m(-0)/m(+0); T -matrix 
with Tn = 1, Ti 2 = and T 2 i ^ corresponds to the 
introduction of a (^-function interface potential; 7 ' 16,17,18 
if only off-diagonal elements contribute to the trans- 
fer matrix we obtain "inverted" boundary conditions 19 
(see also Ref. 5) which hold with a high accuracy for 
GaSb/InAs interface; 5 etc. In general all components 
of the transfer matrix T can be nonzero. They reflect 
internal structure of the heterointerface and cannot be 
expressed in terms of only bulk parameters. For differ- 
ent particular cases components of T-matrix were calcu- 
lated using empirical tight- binding and/or pseudopoten- 
tial approaches. 11,12 ' 13,20 
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The general connection rules Eq. (3) can be rig- 
orously justified within the generalized effective mass 
theory. 21 ' 22,23 ' 24,25 This approach leads to a set of 
integral-differential equations for envelope functions 
which arc defined using a single Bloch basis for the whole 
structure. Coefficients in these equations and, therefore, 
envelopes are smooth and continuous functions even for 
a system with microscopically abrupt interfaces. Near an 
interface the coefficients depend on microscopic details of 
the interface. Normally a perturbation, which is caused 
by the interface, is localized at the atomic scale. Hence, 
if we are interested in behavior of envelop functions on 
the scale which is larger than the lattice constant, we can 
use extrapolated bulk envelopes instead of the exact en- 
velop functions. Despite the exact envelopes are smooth 
and continuous, the extrapolated functions obey general 
connection rules Eq. (3) with parameters which depend 
on details of the interface. 26 The regular calculation of 
the transfer matrix using the generalized effective-mass 
theory 21 ' 22,23 ' 26 is an extremely tedious task. However 
this theory can be considered as a foundation for the 
phenomcnological introduction of the transfer matrix. 

In this paper we follow such a phenomenological ap- 
proach and develop a general method for construction of 
the transfer matrix (Sec. II). Namely, we assume that the 
differential equation (or system of equations) Eq. (1) with 
piecewise smooth coefficients is not applicable at inter- 
face points (which are the points of discontinuity of the 
coefficients in Eq. (1)). The Hamiltonian H in Eq. (1) is 
defined on a space of piecewise smooth and continuous 
functions with linear connection rules at the points of dis- 
continuity. We show that Hermiticity of the Hamiltonian 
on this space of functions impose the first set of restric- 
tions on the form of connection rules. A particular con- 
sequence of these restrictions is the conservation of the 
flow at the interface. Similar methodology, which can be 
found in quantum mechanics text books, 27 has been re- 
cently applied to the one-band effective mass model with 
a general form of the kinetic energy operator. 28 The sec- 
ond set of restrictions follows from the symmetry - the 
transfer matrix must be invariant with respect to trans- 
formations of the interface symmetry group. These two 
sets of restrictions strongly reduce the number of com- 
ponents in the transfer matrix. The rest of T-matrix 
along with band offsets should be considered as empir- 
ical parameters which are defined from experiment (see 
for example Ref. 30) and/or ab initio calculations. The 
method developed in Sec. II is closely related to the com- 
mon method of invariants 29 which allows to construct ef- 
fective k • p Hamiltonians for bulk semiconductors using 
only Hermiticity and symmetry requirements. 

The general connection rules Eq. (3) allow to describe 
various physical consequences of nontrivial internal de- 
tails of a hctcrointcrface. It has been demonstrated in 
Ref. 30 that the use of general boundary conditions re- 
moves quantitative discrepancies between square well cal- 
culations and experiment. 31 There are also more trans- 
parent qualitative effects which come from the complex 



structure of boundary conditions. For example, in multi- 
band models the off-diagonal element T21 (which is equiv- 
alent to an interface 5-potential) is responsible for the in- 
terface heavy-light hole 18 and T — X electron 12,32 ' 33 mix- 
ing. In the present paper we concentrate on a description 
of interface localized and resonant states. 

In 1932 Tamm 34 demonstrated the existence of elec- 
tronic states localized at a surface of a semiconductor. 
It is quite natural to expect that similar localized states 
with energies inside the forbidden gap can occur at an 
abrupt heterojunction. Such a possibility was quali- 
tatively considered by James 35 and later by Zhu and 
Kroemer. 16 In Ref. 36 the existence of interface donor 
states was postulated to explain anomalous transport 
properties of undoped InAs/AlSb quantum- well struc- 
tures. However, general conditions for occurrence of in- 
terface states remained unclear for a long time. This 
problem was addressed in Ref. 37 and recently in Ref. 38 
using the tight-binding approach. Since the effective- 
mass method is of extreme importance for heterostruc- 
ture applications, it is desirable to have a description of 
interface states in terms of envelope functions. 

It is known that envelope-function models with 
the standard boundary conditions possess interface 
states. 39,40 However the corresponding energy levels al- 
ways lie in the region of band offsets outside the energy 
gap of a heterojunction. 39 ' 40 An interesting exception is 
the formation of localized states at heterojunctions with 
band inversion. 39,41 ' 42 These states have a topological na- 
ture and are related to the supersymmetry of an inverse 
contact. 41 

In Sec. Ill of our paper we show that different types of 
interface states which have energies in the forbidden gap 
can be described using generalized connection rules. We 
apply the boundary conditions derived in Sec. II to the 
one-band effective-mass model and to the 4-band Kane 
model which describe T-point states in III-V semicon- 
ductors. We derive general conditions for the existence 
of interface states, and discuss the physical meaning of 
the off-diagonal components in the transfer matrix. In 
Sec. IV we study a scattering problem and demonstrate 
the existence of resonant states and resonant tunneling 
through a single barrier stricture, which are related to 
nonzero off-diagonal elements in the transfer matrix. In 
Sec. V we summarize our results. 



II. GENERAL APPROACH AND BASIC 
EQUATIONS 

To establish a general form of connection rules we con- 
sider the standard statement of a problem within the 
envelope- function approach that is to find eigen functions 
tp(r) and eigen values E of a Hamilton operator H. Let, 
as usual, the Hamiltonian H be a second-order matrix 
differential operator 

H = H Q (r) - ^M a0 d a dp + iL a 8 a (4) 
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where Ho(r), M a p and L a are mxm Hermitian matrixes 
and a, (3 = x,y,z. We assume that the growth direction 
of a structure coincides with z-axis and the system is 
spatially homogeneous in x — y-plane. In this case wave 
functions take the form ip(r) = e ±r ^(z) (kj_ is the mo- 
mentum perpendicular to z-axis). The function ipiz) is 
a solution to the one-dimensional Schrodinger equation 
with the Hamiltonian 



H z = h(z) - 7 9 z 2 + iPd z 



(5) 



where 



h(z) = H {r) - -M a pk± a k±/3 + iL a kj_ c 

7 = M zz , 

P = L z +iM az k ±a 



(6) 



We consider the system which consists of N regions with 
different material parameters and introduce notations 
z n _\ and z n for the left and the right boundaries of the 
nth region (n = 1...N). Assume that h(z) is a piecewise 
smooth function of z, whereas 7 and P are piecewise con- 
stant functions. In the nth region (z n -i < z < z n ) ma- 
trixes h(z), 7 and P respectively take the values h n (z), 
j n and P n , where h n (z) is a smooth and continuous func- 
tion, and 7„ and P n are constants. 

Let Wt be the space of functions ip(z) which arc 
piecewise smooth and square integrable on every inter- 
val z n -i < z < z n . Besides, at every point z n values 
ip(z n + 0), d z tp(z n + 0) and ip(z n - 0), d z ip(z n - 0) arc 
connected by linear relations 



V>(z„ + 0) ^ _ 

d x ll>{Zn + 0) 



^{Zn - 0) 

d z il>(z n - 0) 



(7) 



where T n (n = 1..N) are 2m x 2m matrixes. 

The differential operator H z Eq. (5) is well defined 
on the space Wt, but not necessarily Hermitian. The 
Hcrmiticity condition impose a restriction on a possible 
form of the transfer matrixes T n . 

By the definition the operator H z is Hermitian on the 
space Wt if 



/=< <p\H z \4,> ~ <4,\H Z \^>*=0, 



(8) 



where <p and ip belong to Wt- Matrix elements in Eq. (8) 
are defined as integrals over the whole structure with the 
points of discontinuity being excluded 



< (p\H\ip >- 



N 

E 

n=l. 



(p + (h n (z) - j n d z + iP n d z )ipdz. 



(9) 

We assume for definitcness that zq = —00, z^ = 00 and 
V>(±oo) = 0. It is convenient to introduce 2m-component 
vectors 



d z ip(z) 



<p(z) 
d z tp(z) 



(10) 



and a "current" operator J n which acts on these vectors 



Pn 







(11) 



Integration by parts in Eq. (8) leads to the following ex- 
pression for the quantity / 



N 



7 = z^A„{$+J*}, 



(12) 



71=1 



where A„{$+ J\E'} is a jump of the quantity $+ J* at 
the point z = z n : 



A n {$+J*} 



$+(z„ + 0)J*(z„ + 0) 
$+(z„-0)J*(z„-0). 



The Hermiticity condition I — is fulfilled if 

A„{$+ J*} = 

for every boundary and any pair of functions ip and 
if> from Wt- Using the definition of transfer matrixes 
Eq. (7) we arrive at the Hermiticity condition of the fol- 
lowing form 



Jn T n Jn+lPn? 



(13) 



which means the invariance of the "current" operator J 
under the transfer across a discontinuity point. 

To simplify formulas we consider bellow a system with 
a single boundary at z = which separates left (n = 1 = 
L) and right (n = 2 = R) regions. Hence the connection 
rules Eq. (7) take the form 



tffl(0) = T* L (0), 



(14) 



where the transfer matrix T must satisfy the following 
Hcrmiticity condition 



J r = T + J,T. 



(15) 



Another set of restrictions follows from the fact that 
T-matrix in Eq. (14) should be invariant with respect to 
the symmetry group G of the interface plane: 

D(g)TD- 1 (g) = T, 

where D(g) is 2m x 2m matrix which corresponds to an 
element g of the group G. Since ip(0) and d z tp(0) have 
the same transformation properties with respect to oper- 
ations of the interface symmetry group G, matrixes D{g) 
take a block diagonal form 



D(g) 



D{g) 

D{g) 



(16) 



where mxm matrixes D(g) form a representation (re- 
ducible in general case) of the group G in the basis which 
corresponds to the bulk Hamiltonian Eq. (4) . Hence the 
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FIG. 1: One-band potential profile of the heterojunction (see 
text). 



symmetry requirements can be written independently for 
every mxm block (i, j — 1, 2) of the full transfer ma- 
trix T: 



D(g)T l] D- 1 (g)=T lj 



(17) 



Equations (15) and (17) provide a complete set of phe- 
nomenological requirements which restrict the form of 
the general connection rules Eq. (14). In the next section 
we present a solution of these equations for the one-band 
effective mass model and 4-band Kane model. 



III. LOCALIZED TAMM-LIKE INTERFACE 
STATES WITHIN THE ENVELOPE FUNCTION 
APPROACH 

A. One-band effective mass approximation 

In this section we consider a single heterointerface lo- 
cated at the point z = 0. The corresponding band di- 
agram is shown in Fig. 1. The one-band effective mass 
Hamiltonian takes the form 



V 2 

H n = — h V n , 

2m n 



(18) 



where n = L,R and rriL,R and Vl,r are effective masses 
and potentials in the left (L) and the right (R) bulk re- 
gions respectively. The "current" operator Eq. (11) for 
this model is J n — a y /2m n (a y is the Pauli matrix). 
Hence the Hermiticity condition Eq. (15) takes the form 



I3<i y = T + a v T, p 



m R 
rriL 



The general solution to this equation is 



(19) 



(20) 



where \ is an arbitrary phase and t is a real 2x2 matrix 
with unit determinant: 



dct t = <n<22 - £12*21 = 1- 



(21) 



Let us study localized solutions which are allowed by 
the boundary condition Eq. (14) with the transfer matrix 



Eq. (20). We assume for definiteness that Vl — 0, and 
Vr = V > is a band offset (see Fig. 1). Wave functions 
of interface states take the following general form 



ipR = Ae 



Z>0 
Z<0. 



ipL = Be ik±r e KLZ , z<0. (22) 
The energy of the localized state is defined as 



E 



2uil 2tol 



where K 2 L /2mL is the binding energy. Quantity kr in 
Eq. (22) is related to kl by the equation 



K R = ^pn 2 L +q 2 , q 2 = 2m R V+(l-l3)k 1 _ 2 . (23) 

Parameter q in Eq. (23) describes a degree of the het- 
erointerface asymmetry which is related both to the band 
offset and to the difference of effective masses. 

Substitution of Eq. (22) into the boundary conditions 
leads to the following dispersion equation for kl 



(3*1 + q 2 



*21 + *22^L 
til + ti 2 K L 



(24) 



It is natural to assume that the diagonal elements of 
the matrix t are positive (in > 0, £22 > 0). In this 
case Eq. (24) has real solutions only if £21 and/or tvz are 
negative. There are two types of solutions. Solutions of 
the first type correspond to the case £21 < 0, t\i > 0. 
The existence condition takes the form 



*21 

— >q- 

Ell 



(25) 



These states occur near the lowest band edge and move 
down into the energy gap with increase of |*2i|- Equa- 
tion (25) shows that the existence of such solutions is 
restricted by the value of the asymmetry parameter q 
Eq. (23). It was mentioned in the Introduction that a 
nonzero element £21 can be modeled by an interface <5- 
potential. Negative £21 corresponds to an attractive in- 
terface potential. Hence Eq. (25) is analogous to the well 
known condition for the existence of bound states in a 
potential well with asymmetric barriers. 43 If ^ tur 
{[3 7^ 1) the asymmetry parameter q depends on 2D mo- 
mentum k± (see Eq. (23)). Therefore Eq. (25) defines a 
line in k^-space which separates localized 2D interface 
states and delocalized 3D continuum states. Analogous 
2D-3D transformations were recently studied in asym- 
metric quantum wells 44 . 

Solutions of the second type are related to negative 
values of the second off-diagonal element ti2 . They exist 
at arbitrary £12 < and £21 > 0. At small negative £12 the 
energy levels, which correspond to these states, lie deep 
in the forbidden gap and approach the band edge with 
increase of |*xa ] - The solutions of the second type can be 
naturally viewed as states which originate from a lower 
(for example valence) band and move up with increase 
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of an effective interface potential. In Sec.IIIC we shall 
return to this point and discuss a possible interpretation 
of the element ti2 in terms of a local perturbation of 
remote bands. 

Obviously, the solutions of both types can coexist if 
both off-diagonal elements are negative. 



B. Interface states in the 4-band Kane model 

To study interface states in multi-band systems and to 
illustrate the importance of the symmetry requirements 
Eq. (17) we consider the 4-band Kane model which de- 
scribes T-point states in III-V zinc blend semiconductors 
without spin-orbit splitting. We consider a single (001) 
heterojunction and find the interface states which origi- 
nate from A-line states with kj^ = 0. 

It is convenient to introduce the following basis 



{ Xj } = {\S),\Z),\X + ),\X_)}, 



(26) 



where \X±) = -±=(\X) ± \Y)). In this basis the A-line 
Hamiltonian takes a block diagonal form 29 



HI 



(27) 



where H" z and H± are 2x2 matrixes which correspond 
to \S), \Z) and \X + ), pairs of states respectively 

(n = L, R) 



Tjn 



HI = 



P d 

2m- r n<~> z 



E vn ~\~ 



E - 



2m v 



(28) 



In Eq. (28) E cn and E vn are the energies of conduction 
and valence band edges respectively and / is the 2x2 
unit matrix. 

First we establish the general form of boundary condi- 
tions for this system. The transfer matrix 



T = 



Tn T12 
T21 T22 



(29) 



which enters the connection rules Eq. (14) consists of four 
4x4 blocks. Every block Tjj = 1,2) must satisfy 
the symmetry conditions Eq. (17). The symmetry group 
of (001) plane for zinc blend structure is the group C^ 18 
which has four elements: the unit element E, a second- 
order axis C2 and two mutually perpendicular reflection 
planes 01 and 02- This group has four classes and thus 
four irreducible representations. Each function from the 
set {Xj} Eq. (26) is the basis function for one of the ir- 
reducible representations. Namely, functions \S) and \Z) 
correspond to the same representation A x , whereas func- 
tions \X + ) and \X_) are related to the representations 
Bi and B2 respectively. Thus matrixes D(g) which enter 



the symmetry conditions Eq. (17) have a diagonal form 
in the basis Eq. (26) 



D(E) = diag(l, 1,1,1), 

D(C 2 ) = diag(l, 1,-1,-1), 

D(<7i) = diag(l, 1,1,-1), 

D{a 2 ) = diag(l, 1,-1,1), 



(30) 



where diag(...) stands for the set of diagonal elements. 
Straightforward calculations show that the general solu- 
tion to Eq. (17) with D(g) Eq. (31) takes the form 



TH 
T| j 
T ij 



(31) 



where T l J z is an arbitrary 2x2 matrix and T± are ar- 
bitrary numbers. Consequently, pairs of states (\S), \Z)) 
and (|X + ), |X_)) as well as the states \X + ) and \X_) 
are decoupled due to the symmetry requirements. Thus, 
we get three independent problems. Two of them cor- 
respond to the solution of two independent one-band 
Schrodingcr equations with general boundary conditions 
(see Sec. IIIA) for decoupled \X + ) and \XJ) valence 
bands. It is worth to mention that if TV ^ Tj 7 , which is 
allowed by the symmetry, the connection rules Eqs. (14), 
(29), (31) lead to a heavy-light hole mixing. In fact, the 
connection rules 



4>x,yl(0) = *Px,yr(0), 



(32) 



which arc used 18 to describe the heavy-light hole mixing 
at the normal hole incidence, represent a particular case 
of the T-matrix Eqs. (29), (31) with T^ = 1, T^ 2 = 
and T 21 = -T 21 = T^ h . 

To describe localized states which correspond to the 
subspace {|5'),|Z)} we have to solve the two-band 
Schrodinger equation with the Hamiltonian H sz Eq. (28). 
The boundary conditions to this equation are defined via 
the transfer matrix T sz which has no symmetry restric- 
tions since both \S) and \Z) correspond to the same rep- 
resentation Ai of the interface group C2V 

To simplify further calculations we neglect the second- 
derivative terms in H sz . This reduces the problem to the 
solution of two coupled first-order differential equations 



E cn — E P n d z 

Pn&z E vn E 



V>„(z) = 0. 



(33) 



Since the highest spatial derivative in Eq. (33) is of the 
first order we should not include the derivatives of wave 
functions in the boundary condition. Therefore the trans- 
fer matrix T sz has only one nonzero 2x2 block T^ z = T sz 
which is restricted only by the Hermiticity condition. 
The "current" operator for the problem Eq. (33) takes 
the form Jj = Pjcr y . Hence the Hermiticity condition 
Eq. (15) formally coincides with Eq. (19) 



7 = 



Pl 
Pr 



(34) 
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(a) 




FIG. 2: (a) Band digram for a hetero junction of the I type; 
(b) Band digram for a heterojunction of the II type 



The solution to Eq. (34) is an arbitrary real 2x2 matrix 
with fixed determinant and with an arbitrary phase (see 
Sec. IIIA) 



T sz =te lx ^j, det< = l. 
Considering a localized interface solution 



(35) 



UL 
VL 



we get the following dispersion equation 

ti2-tnF L (E) 



Fr(E) = 



to 



t 2 iF L {Ey 



where 



F n (E) = 



E — E„ 



En 



E 



(36) 
(37) 

(38) 



The localized solution of the form Eq. (36) exists if the 
energy E lies in the forbidden gap of the heterojunction 

max{E V L, E v r} < E < min{E C L, E c r} 

which is shown by shaded region in Fig. 2. It is conve- 
nient to introduce a new variable 



x{E) = F L {E) 



E-E 



vL 



(39) 



E c l — E 

and rewrite the dispersion equation Eq. (37) in the form 

ti2 - t\\X 



Fr(x) 



t 22 - t 2 lX 



(40) 



where the function F R {x) is defined as follows 



P / \ E vL -E vR + x 2 (E cL -E vR ) 

y E cR — &vL + X^yticR ~ &cL) 

We analyze the solutions of Eq. (40) under the natural 
assumption t\x > and £22 > 0. We also assume that 
E c r > E c l i.e. the conduction band offset V = E c r—E c l 
is positive (see Fig 2). Therefore a heterointerface of the 
I type (Fig. la) corresponds to the condition 



E v r < E v l, 
whereas the inverse inequality 

E v r > E v l 



(42) 



(43) 



holds if the heterostructure belongs to the II type (see 
Fig. 2b). 

First we consider heterointerfaces of the I type. In 
this case the region of the energy gap (shaded region in 
Fig. 2a) 



E v l < E < E c l 



maps to the region 



< x < 00 



of the variable x. When x goes from to 00, the function 
Fr(x) in the right hand side of Eq. (40) monotonically 
increases from 



^nin = V(E vL ~ E vR )/(E cR - E vL ) (44) 
at x = (which corresponds to E — E v l) to the value 

FLx = V(E cL - E vR )/(E cR - E ct ) (45) 

at x = 00 (E = E c l). The function Fr(x) for a hetero- 
junction of the I type is shown by solid line in Fig. 3. The 
behavior of the right hand side in Eq. (40) depends on the 
signs of the off-diagonal elements £12 and t 2 i- The dis- 
persion equation Eq. (40) have real solutions if at least 
one of the off-diagonal elements is positive. If tyi and 
£21 have opposite signs there exists only one solution to 
Eq. (40). 

Let us analyse different cases separately. 

(i) £12 > 0, £21 < 0. In this case the right hand side in 
Eq. (40) is a decreasing function of x. A solution exists 
if 



h 

t 22 



> 



E v l — E v r _ 



E, 



cR 



E vl 



(46) 



This solution can be identified with an acceptor interface 
state which originates from the valence band and move 
up from the edge of the valence band E v l to the edge of 
the conduction band E c l with increase of t\2- 
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FIG. 3: Dependence Fr vs x (see Eq. (40) in the text). Solid 
(dotted) line shows Fr(x) for the heterojunction of the I (II) 
type. 



(ii) i 2 i > 0, t 12 < 0. A solution to Eq. (40) exists 
under the following condition 



hi 
tn 



> 



E, 



cR 



E, 



cL 



E r 



E. 



vR 



= (f i y 1 

V max/ 



(47) 



This is a donor interface state originated from the con- 
duction band. With increase of t 2 i the energy level, 
which corresponds to this state, moves through the en- 
ergy gap from E cL to E vL . 

If the heterointerface belongs to the II type, then E v r 
forms a lower bound of the energy gap (shaded region in 
Fig. 2b). Hence the energy gap 

E v r < E < E cL 

maps to the region 

x min = V (E vR - E vL )/ (E cL - E vR ) < X < oo 
in the x-axis. 

The function Fr(x) in Eq. (40) increases from zero at 
x = x min (E = E vR ) to the value i^/ ax = F^ ax Eq. (45) 
at x — oo (E — E c x). This function is shown by the 
dotted line in Fig. 3. Consequently, only the condition 
for the existence of the acceptor state is changed. For 
the heterointerface of the II type the existence condition 
Eq. (46) is replaced by the inequality 



*ia 
tn 



> 



E v r — E v l 



E, 



cL 



E V R 



(48) 



The condition for the existence of the donor state Eq. (47) 
remains unchanged. 

Acceptor and donor interface states coexist if t 12 t 2 i < 
iii*22 (*12 > 0, far > 0) and the conditions Eq. (46) 
(Eq. (48) for the II- type structure) and Eq. (47) are ful- 
filled simultaneously. 



Physical meaning of the off-diagonal elements in 
the transfer matrix 



In this subsection we discuss a possible interpretation 
of the elements T\ 2 and T 2 \ in the transfer matrix for the 
one-band effective mass model. 

It is well known 17,20 that the element T 2 \ can be in- 
terpreted as an interface (^-function potential. Indeed, 
a nonzero T 2 \ introduces a jump of the first derivative 
which is proportional to the value of the wave function, 
exactly as an interface (^-potential does. The physical 
meaning of the second element T 12 is less clear (see, for 
example, discussions in Refs. 17,20,26). In Sec. IIIA we 
have shown that there exist localized states which are re- 
lated to nonzero T\ 2 . These states behave as "acceptor" 
states which originate from remote lower bands. There- 
fore, it is natural to expect that the element T\ 2 is related 
to a local perturbation of these remote bands though they 
are not explicitly included in the bulk one-band Hamil- 
tonian. To confirm this interpretation we consider a one- 
band model as the limit of the two-band model with the 
following Hamiltonian 



H 



E c (z) Pd z ' 




' 9c " 


~P3 Z E v (z) _ 


+ 


g v 



6(z), (49) 



where E c t v \ (z) describes the profile of the conduction (va- 
lence) band. The first term in Eq. (49) corresponds to 
bulk regions whereas the second term models a local per- 
turbation which is caused by the interface. 

The boundary conditions for the two-component enve- 
lope function can be obtained by the integration of the 
Schrodinger equation 



H 



= E 



over a small segment [—a, a] and taking the limit a 
The result of the integration takes the form 



(50) 
-> 0. 



M0) = §Iim 

f- a— *0 



Mo)-^(o) 

Using the identity 



— lim 

P a^O 



8{z)6{z)dz 



8{z)(p{z)dz, 

a 

6(z)ip(z)dz. 



(9{z) is the Heaviside function) we arrive at the connec- 
tion rules 



where 



Vr 



tu — til — 



t 2 \ t 22 



1 - g c g v /AP 2 
1 + g c g v /APi 



9v/P 



9c/P 



l + g c g v /AP^ 



1 + g c g v /4P- 



(51) 

(52) 
, (53) 
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Equations (51)-(53) show that both off-diagonal elements 
for two-band model are reproduced by the simple inter- 
face term Eq. (49), though the transfer matrix in Eq. (51) 
is still not of the most general form (compare to Eq. (35)). 
We can also clarify the physical meaning of the interface 
solutions for the two-band model which has been consid- 
ered in Sec. IIIB. According to the results of Sec. IIIB the 
interface solution of acceptor type exists if ti 2 is positive. 
Positive ti2 corresponds to positive g v (see Eq. (53)) and, 
consequently, to a local perturbation of the valence band 
which is attractive for holes. Analogously, the solution of 
the donor type is related to a positive value of £21 which 
corresponds to negative g c and a local perturbation of 
the conduction band, attractive for electrons. 

Let us derive the one-band model which is related to 
the two-band Hamiltonian Eq. (49). If the energy E in 
Eq. (50) is close to the edge of the conduction band 

\E-E c (z)\/A(z)^l 



IV. RESONANT TUNNELING THROUGH A 
SINGLE BARRIER WITH COMPLEX 
INTERFACES 

As we have seen in Sec. Ill, a nontrivial internal struc- 
ture of a single heterointerface which is described by the 
generalized boundary conditions allows for the existence 
of localized interface states. In this section we show that 
the interference effects in a system with more than one 
interface can result in the formation of resonant interface 
state. The most transparent manifestation of these states 
is the resonant tunneling through a single-barrier struc- 
ture. We shall demonstrate this effect for the one-band 
effective mass model. 

Let us consider a single-barrier heterostructure with 
a rectangular potential barrier of the height V and 
the width L. The standard scattering solution to the 
Schrodinger equation is defined by the following asymp- 
totic form of the wave function 



(2 A = E c — E v is the energy gap), we can express the 
lower component of the spinor in Eq. (50) in terms of the 
upper component 



<p(z) 



2A(z) 



d z tp(z). 



(54) 



The upper component ip(z) plays the role of the wave 
function which in the bulk regions satisfies the one-band 
Schrodinger equation 



(E — E c (z))ip(z), 



(55) 



where m(z) — A(z)/P 2 is the effective mass. Equa- 
tion (55) should be supplemented by boundary conditions 
which are obtained by the substitution Eq. (54) to the 
connection rules Eq. (51). The final boundary conditions 
for the one-band model take the form 



Tn T12 \ / ip L 
T21 T 22 ) V 9 z tp L 



(56) 



with the following elements of the transfer matrix 



rp _m L _ 1 - 9c9v/^P 2 
m R l+g c g v /4:P 2 



T12 = - 



1 + g c g v /AP 2 ' 



T21 = 



^m R g c 



l+g c g v /AP 2 



(57) 
.(58) 



Thus, the off-diagonal element T12 Eq. (58) is propor- 
tional to the strength g v of the local perturbation of the 
remote valence band. A potential, which is attractive for 
holes, corresponds to positive g v and thus negative Ti 2 . 
This explains the results of Sec. IIIA and confirms our in- 
terpretation of the interface state related to the element 
T12. 



rue 



-ikz 



t k e 



ikz 



— OO 
+OO 



(59) 



where the wave vector k is related to the energy of the 
incident wave E = k 2 /2m. Using the general connec- 
tion rules Eq. (14) and performing straightforward cal- 
culations we arrive at the following expression for the 
transparency coefficient Dk — \tk\ 2 

D k = Ak 2 K 2 { [(k 2 T^ + k 2 Tl 2 + T 2 2 ! + k 2 n 2 Tl 2 ) sinh kL 

+ 2K(T 11 r 21 + fc 2 T 22 r 12 )cosh K L] 2 + 4fc 2 K 2 } 1 (60) 

where 



k = \j2mV - k 2 . 

The transparency Dk Eq. (60) has an explicit resonant 
structure and turns into unity at the condition 



tanh K L= — ^ TllT21 + fc2722Tl2 ) 

k 2 ? 12 ! + k 2 T 2 2 + T 2 i + k 2 n 2 T 2 2 



(61) 



If L — > 00, Eq. (61) has no solutions. However, the reso- 
nance occurs if the width of the barrier L becomes smaller 
than some critical value L c and if at least one of the off- 
diagonal elements is negative. 

To reveal the physical nature of the resonance condi- 
tion Eq. (61) we consider the simplest nontrivial trans- 
fer matrix with only one nonzero off-diagonal clement. 
Namely, we assume that Tn = T 22 = 1, Ti 2 = and 
T 2 i 7^ 0. This transfer matrix corresponds to the (5- 
function interface potential of the strength g = 2toT 2 i. 
Since the resonant solution exists only for negative T 2i we 
also assume that T 2 i < 0. Under the above assumptions 
the resonance condition Eq. (61) reduces to the equation 



tanh kL 



2k\T, 



21 



2mU 



t 2 Y 



(62) 



FIG. 4: Functional dependence of r on \T2i\/V2mV 
(Eq. (64)). 







1 



x 



Introducing the following dimcnsionless variables 



FIG. 5: Graphical solution of the equation tanhxl = tx 
(Eq. (66)) which defines the behavior of the under-the-barrier 



K I E 

„ \T 21 \/V2mV 



1+ (\T 21 \/V2mV 



I 



LV2mV, 



we transform Eq. (62) to the form 
tanhxl = tx. 



(63) 
(64) 

(65) 
(66) 



Parameter t in the right hand side in Eq. (66) can not ex- 
ceed unity (r < 1). In fact, the function r(|T 2 i \/V2mV) 
Eq. (64) (see Fig. 4) reaches its maximum value r max = 1 
at 



\T 21 \/y/2mV = 1. 



(67) 



Note, that Eq. (67) defines the critical value of IT21 1 
which is required for the existence of the localized state 
at the single heterointerface. Indeed, Eq. (25) shows 
that the interface bound state exists if \T 2 i\/\ / 2mV > 1. 
Since the behavior of solutions to Eq. (66) is governed by 
the parameter r, the existence of the under-the-barrier 
resonances does not depend on whether or not the sep- 
arate single interface possesses bound states. We shall 
see, however, that these resonances are closely related 
to the bound states which exist in the combined double- 
interface structure. 

Under-the-barrier resonances correspond to the solu- 
tions to Eq. (66) in the region 0<x<l(0<E<V). 
Graphical solution of Eq. (66) is illustrated in Fig. 5 
where the functions tx and tanh xl are plotted by dashed 
and solid lines respectively. The required solution exists 
if parameter I, which is proportional to the width of the 
barrier, satisfies the inequalities 

t <l <l c , (68) 

where l c is the solution to the equation tanh l c = t 

l c = -ln^l. (69) 
2 1-t v ' 



resonance. 



Figure 5 shows that with decrease of I we first meet the 
resonance condition at I = l c . The solution x = 1 cor- 
responds to the resonant transparency Dk = 1 at zero 
energy E = 0. With the further decrease of I the reso- 
nance energy moves up, reaches the the top of the barrier 
E = V (x = 0) at I = t, and enters the over-the-barrier 
continuum. 

The asymptotic form of the scattering solution Eq. (59) 
at k = — iq, ti q — ±1 and Ti q = exactly coincides 
with the asymptotic form of a localized solution. Hence 
the condition for the resonant transparency at zero en- 
ergy |ife=o| = 1 is, in fact, the condition for appear- 
ance/disappearance of a localized state (symmetric or an- 
tisymmetric). Thus, the resonant state is nothing but the 
localized state (obviously antisymmetric) which is pushed 
out of the band gap. 

Now we are able to describe the evolution of both res- 
onant and localized states with the change of the barrier 
width L at fixed V and |T 2 i|. We consider separately two 
different cases: 

(i) \T 2 i\/V2mV > 1. At L — > oo there exist two de- 
generate localized states which arc related to the well 
separated heterointerfaces. With decre ase of L the de- 
generacy is lifted and, at L = L c = l c /\j2mV , the upper 
state is pushed out of the band gap to form the resonance. 
With the further decrease of L the resonance moves up 
and crosses the top of the barrier at L = Ty / 2mV. At 
L < T\f2mV we have no resonance, but only the local- 
ized state. 

(ii) \T 2 \\j\[2mV < 1. At L > L c there exists nei- 
ther localized nor resonant state. The critical value 
L = L c corresponds to the creation of the resonance- 
bound state pair. When L becomes smaller than L c the 
resonance moves up and reaches the top of the barrier 
at L = T\ / 2mV, whereas the localized state moves down 
into the energy gap. 

Thus in either case the resonant transparency is always 
accompanied by the bound state. The resonance and the 
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localized state can be qualitatively interpreted as anti- 
bonding and bonding orbitals respectively. 

It is interesting to note that in the simple case of 
zero element 7i 2 the resonant transparency of the single- 
barrier structure has a clear counterpart in the classical 
electrodynamics. The corresponding system is a metallic 
slab which is covered on either side by dielectric layers. 
The metallic region with e(uj) < models the barrier, 
whereas the dielectric coating corresponds to the attrac- 
tive interface potential. The resonant transparency of 
electromagnetic waves through such a system has been 
described theoretically 45,46 and observed experimentally 
in Ref. 46. If both off-diagonal elements are nonzero, 
the resononant tunneling through a single-barrier het- 
erostructure apparently has no an optical analog. 

V. CONCLUSION 

It it commonly accepted that effects of a microscopic 
structure of a heterointerface can be incorporated into 
the envelope function method by the use of generalized 
connection rules. Such connection rules are normally for- 
mulated in terms of the interface transfer matrix. In this 
paper we presented the general method which allows to 
construct the transfer matrix for arbitrary system. We 
showed that the Hermiticity of the Hamiltonian and the 
symmetry of the interface plane impose the restrictions 
on the form of transfer matrix. These restrictions can 



be formulated as a set of equations for components of 
the transfer matrix. Solution of these equations defines 
the general form of the boundary conditions for a given 
interface. To illustrate this approach we considered the 
one-band and 4-band envelope function models and es- 
tablished the conditions for the existence of Tamm-likc 
localized interface states. We have also demonstrated 
that in a system with more than one heterointerface there 
exists a possibility for new physical effects such as the 
resonant tunneling through a single potential barrier. 

In this paper we did not discuss the application of 
our general results to particular heterostructures. Such 
an application requires a knowledge of phenomenological 
parameters which enter the generalized boundary condi- 
tions. Identification of these parameters using the results 
of experiments or ab initio calculations is by far not sim- 
ple task. We believe, however, that the physical effects 
which have been described in our paper and which origi- 
nate solely from the nontrivial structure of the boundary 
conditions, could be helpful in resolving this problem. 
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